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^ ■ Abstract 

> : 

^ ■ It is suggested that the measurements of hadronic invariant mass (nix) distributees 
in the inclusive B — > X c ^lv decays can be useful in extracting the CKM matrix ele- 
ment \V u b\. We investigated hadronic invariant mass distributions within the various 
theoretical models of HQET, FAC and chiral lagrangian as well as ACCMM model. 
It is also emphasized that the mx distribution even at the region mx > mo in the 
inclusive b — > u are effetive in selecting the events, experimentally viable at the future 
asymmetric B factories, with better theoretical understandings. 
PACS numbers: ll.10.Jj, 11.90.+t 



1. Introduction 



The CKM matrix element V u b is important to the standard model description of 
CP violation. If it were zero, there would be no CP violation from the CKM matrix 
element (i.e. in the standard model), and we have to seek for other source of CP 
violation in Kl — > vr7r. Observations of semileptonic b —>■ u transitions by the CLEO 
U and ARGUS |§ imply that V u b is indeed nonzero, and it is important to extract 
the modulus \V u b\ from semileptonic decays of B mesons as accurately as possible. 
Presently, the charged lepton energy is measured, and the b — ► u events are selected 
from the high end of the charged lepton energy spectrum. This method is applied to 
both inclusive and exclusive semileptonic B decays. 

However, this cut on Ei is not very effective, since only below 10% of b — > u events 
survive this cut. As first discussed in ||, the measurements of hadronic invariant mass 
(mx) distributions in B — > X c ^ u lv (inclusive decays) can be useful to extract a CKM 
matrix element V u b- For b — > civ, one necessarily has mx > % = 1.86 GeV. So, if we 
impose a condition mx < mn, the resulting events come from b — > ulv. According to 
the work of [[J, ~ 90% of the b — > u events survive this cut. This is in sharp contrast 
with the usual cut on E\. Thus, one could get an order of magnitude improvement in 
selecting data. 

Inclusive mx distributions for b — > u transition can be obtained using the ACCMM 
model 0, or the work by BUSV ||. In this work, we follow the result of Ref. || 
based on the ACCMM model, with a remark that theoretical uncertainties for the 
inclusive b — > u decays are less compared to an exclusive mode. 

Resonance contributions to the mx distributions in B — > X CjU lv can be estimated 
invoking various models. Once the decay rate for B —> Rlv (where R is a resonance) 
is known, the corresponding mx distribution can be written as 

dT ^ 2m x T(B -» Rlv) m R T R 

dmx 7T ( m x ~ m i?) 2 + m R^ 2 R 

in the narrow width approximation. Here, m R and Tr are the mass and the width of 
the resonance R. In the limit of Tr — > 0, we get 

= T(B -> R) 5(m x - m R ), (2) 

dmx 

so that the correct decay rate for B — > Rlv comes out upon integrating Eq. (2) over 
dmx- 
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In Section 2, we discuss the result for B — > D,D*,(D**) in the heavy quark 
effective theory. In Section 3, B — > (jr,p)lu are considered in two different types 
of approaches: a nonrelativistic quark model and the chiral lagrangian with heavy 
mesons as well as light vector mesons. In Section 4, the mx distributions for B — > 
X c < u \lv are shown, and it is emphasized that the mx distribution even at the region 
mx > tud in the inclusive h — > u are effetive in selecting almost 100% of the events, 
experimentally viable at the future asymmetric B factories. Since one can calculate 
the inclusive decay more reliably, one can achieve better determination of V n b both 
statistically and systematically. 



2. The B —> X c lv decay in the heavy quark effective 
theory 

Let us first consider B — > X c lu, which is dominated by resonance contributions 
with X c = D, D*, D**. Theoretical predictions based on the heavy quark effective field 
theory (HQET) || depend on one hadronic form factor h^iw). In order to calculate 
a decay rate, one has to know the shape of the form factor over the whole kinematic 
range of w. However, this form factor is not calculable from the first principle, except 
for the normalization : ^Ai(l) = 0.99 ±0.04, which is one of the predictions of HQET 
f|. Thus, one necessarily resort to some model for the shape of the form factor. If 
one adopts the result of the QCD sum rule results, one may approximate the form 
factor as 

h Al (w)a 0.99 (— — ) , (3) 

with q\ pa 0.8. Since this result is based on the QCD sum rule, there exist some 
systematic uncertainties associated with the sum rule. This systematic uncertainty 
may be taken into account by allowing q\ vary between 0.5 and 1.1, where the latter 
is on the border of the limit given by Voloshin's sum rule 0. For this range of g\ , 
one can predict the decay rates for B — > D^*'lv : 

F(B -> Dlv) = (0.86 ~ 1.35) x 10 13 |H c | 2 /sec, (4) 
T(B -> D*lv) = (2.59 ~ 3.48) x 10 13 |V fec | 2 /sec, (5) 

where the smaller decay rates correspond to the larger g\ . For B — ► D**lu, we use 
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the observation by CLEO ||: 

Br(B -> £>**///) ~ 0.5 Br(B -»• DZi/). (6) 

Our approach concerning the measurement of |V^| from the distributions can 
be regarded independently of the uncertainties in Eqs. (4)-(6), because the m x 
distributions from the three resonances D, D* and/?** are essentially delta functions, 
Eq. (2), and we just require to exclude small regions around mx = mp, m^., m£>**. 
For more details, see Section 4 and Fig. 1. 



3. The B X u lu decay 

Unlike the b —>■ c transition, the b —>■ u transition is largely nonresonant and mul- 
tiple jet-like final states dominate 0. The whole inclusive decay can be theoretically 
well understood in most of the kiematical region. The electron energy spectrum or 
the hadronic mass distribution for the inclusive semileptonic decay can be calculated 
rather reliably. In contrast, for the exclusive decays for b — > u such as B — > (jr, p) + lv, 
the model dependence becomes more pronounced, especially for the shape of the form 
factors. Here, we consider two classes of models, the FAC model (a nonrelativistic 
quark model) and the chiral lagrangian with heavy mesons. The results are compared 
with the mx distribution obtained by the ACCMM model in Section 4. 



A. The FAC model 



The FAC model is based on the nonrelativistic quark model with assuming the 



form factors are factorized as |TH] 



Zf A V) = /f Q V) x F(q% (7) 

where / J F ^ M ((? 2 ) is the free quark model form factors arising from the overlap of the 
spin wavefunctions, and F(q 2 ) comes from the overlap of the spatial wave functions. 
For B — » transitions, one can impose the normalization condition, F(q^ nax ) = 1, 
by heavy quark flavor symmetry. For B — > 7r(or p) transitions, this normalization is 
not valid and it may be reasonable to assume that the form factor suppression for 
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B — > 7r begins at the B ^ p threshold. With these assumptions, one gets 



Y{B - 


-> DZi/) 


= (0.71 ~ 


0.85) 
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10 13 


\V bc 2 /sec 
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2.44) 
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10 13 


\V bc \ 2 /sec 


T(B° - 
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10 13 
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= (0.77- 


2.10) 
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10 13 


\V ub 2 /sec 



(9) 
(10) 

(11) 



for certain ranges of pole masses (see Ref. [I0~| for more details.) Note that the FAC 
model predictions for B — > D* are consistent with (although they are systematically 
lower than) those by the HQET discussed in the previous section. Hence, the FAC 
model for the B — > transitions are rather reliable. For B — > 7r(or p) transitions, 
the predictions are very sensitive to the shape of the form factors because of the 
large phase space available. Therefore, we simply regard the above numbers for 
B — > 7r(or p) transitions as exemplary values in a nonrelativistic quark model, without 
giving much meanings to the specific values. 

B. The chiral lagrangian with heavy mesons 

Recently, the chiral lagrangian with heavy mesons and baryons has been developed 
||11||-||13||. This lagrangian was originally invented in order to describes interactions 
among heavy mesons and light mesons such as tt and K in the soft pion limit. Then, 
heavy baryons [14]] as well as p [ 15 1 — [ 16 1 have been incorporated in the leading order 



in 1/m.Q and chiral expansions. The weak current can be represented in terms of 
physical fields like heavy hadrons and light mesons, allowing us to calculate the matrix 
element of the weak current between hadrons and thus the semileptonic decays of 
heavy hadrons. 

However, this apporach has its own limitations. First of all, the hadronic form 
factors given by this chiral lagrangian is valid only in very limited regions of the 
whole kinematic region. Therefore, one often assumes certain shape of form factors 
and normalize them at a point to a value given by the chiral lagrangian with heavy 
hadrons. Furthermore, if one considers the next-to-leading order corrections in I/tuq 
and chiral expansions, there come in a lot of unknown parameters and one essentially 
loose predictability. Although the reparametrization invariance of the heavy quark 
field leads to some constraints to the parameters in the next-to-leading order terms, 
it still leaves many other parameters unconstrained. Therefore, results based on the 
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chiral lagrangian with heavy baryons should be understood, keeping in mind the 
uncertainties just mentioned above. 

One of the extensive studies of semileptonic decays of heavy mesons in the frame- 
work of the chiral lagrangian with heavy hadrons is the work by R. Casalbuoni and 
his collaborators [15]. Their results are 

2 



r(5 c 



7T 



-lu) 



38.8 



7fl(MeV) V 
. 200 



x 10 13 |Kfe|7sec, 



T(B°^p-lu) = 22.7 (^|^) 2 x 10 13 |K fe | 2 /sec 



(12) 
(13) 



At this point, a remark on the B— meson decay constant in Eqs. (|T2|) and ([13]) is 
in order. In the lowest order in the 1/tuq expansion, 
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Id 



<M 



D 



(14) 



On the other hand, the lattice QCD and the QCD sum rule |L7|] suggest that 

f B « f D « 200 MeV, (15) 



are ex- 



which violates the scaling relation, Eq. (0). Thus, the results in Ref. [15 
pressed as above, although it is not systematic in I/tuq expansion to use Eq. (p~5|) . 

We note that the results of Ref. [15] are subtantially larger than those based on 
the FAC model. Especially, relative ratios between B — > tt and B — > p are opposite 
in two models, and may be checked in the near future. For the isospin-related decay 
B~ — ► p l~vi, the predicted decay rate is the half of Eq. flilf ), with the corresponding 
branching ratio 
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Br(B~ -> p°rvi) = 0.44 x 10" 



MMeV) 
200 



V, 



ub 



0.0045 



(16) 



assuming t b = 1.29 ps. The data from ARGUS and CLEO seem contradictory with 
each other : 



Br(B~ — > p°l~v{) = (1.13 ± 0.36 ± 0.26) x 10 
< 0.3 x 10~ 3 
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(ARGUS) 
(CLEO). 



Note that two data are incompatible with each other. The ARGUS result [18 



;i7) 

'18) 
is 



consistent with the prediction by R. Casalbuoni et ai, but is inconsistent with the 
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FAC model prediction. On the other hand, if the result by CLEO [|TJJ is confirmed, 
the prediction based on the chiral lagrangian with heavy mesons would be excluded. 
In this case, there can be many possible reasons for it. First of all, interactions 
between p and heavy mesons may not be well described by the chiral lagrangian in 
the lowest order because of relatively heaviness of p. This would be contrary to the 
better known case, the chiral lagrangian with vector mesons (p), where dynamics of 
7r, p are rather well described. Secondly, and most likely, the usual simple assumption 



on the shape of the form factor may not be right. In most cases including Ref. [[15 
it is assumed that a form factor f(q 2 ) satisfies a monopole form : 

where M po i e is a pole mass. This extrapolation of the form factor through the whole 
kinematic range does not have justifications from the first principle, and is a source 
of uncertainties in any models. 



4. Discussions and conslusions 

The resulting mx distributions for B — ► Rlv for R = tt, p, D, D* , D** are shown 
in Fig. 1, along with the inclusive mx distribution for the b — > u transition, with 
\Vub/Vcb\ = 1- The b — > c transition is dominated by the X c = D,D*,D**, and can 
be reliably calculated in the HQET as described in the previous section. The regions 
between the triangles are the range of the predicted rate when the dmx integration 
over the delta function is performed. On the other hand, the b — > u transition is 
largely nonresonant. The cases with X u = ir,p are shown explicitly for two different 
models discussed in Section 3. For X u = n, the region between the open triangles are 
predictions by Hagiwara et al. TO], and the region between the closed triangles are 



predictions by Casalbuoni et al. [|L5]. For X u = p, the regions between the lower two 
curves are predictions by Hagiwara et al, whereas the region between the upper two 
curves are predictions by Casalbuoni et al. The inclusive mx distribution for b — > u 
was obtained from the ACCMM model with hadronic mass constraint of mx ~ 2m w . 
The excusive decay of B — > tcIu is shown separately. 

From Fig. 1, most of the b — > u transition events survive the cut on the hadronic 
invariant mass, mx < contrary to the more conventional cut on the electron 
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energy. In fact, one can relax the condition mx < m^, because the nix distribution 
in b — > civ is completely dominated by contributions by three resonances D, D* and 
D**, which are essentially like delta functions, Eq. (2). In other words, one can use 
the b —> u events in the region even above mx = tud, excluding small regions in 
nix around nix — nip, nip*, nip**. The cut on the nix is more effective than the 
cut on the electron energy by factor of ~ 10, and therefore we have much better 
statistics. Furthermore, theoretical understanding of exclusive decay modes of B — > 
X u lv is rather poor, as we discussed in Section 3. Two different models lead to vastly 
different predictions for X u = it and p. This would induce theoretical uncertainties 
in determination of V u b from the measurement of an exclusive semileptonic decay of 
B mesons. On the other hand, the inclusive decay is better understood, so it would 
be more reliable to calculate the nix distribution for inclusive b — > u transitions. 

At future B— factory with microvertex detectors (symmetrical or asymmetrical), 
one expects that the efficiency for the event reconstruction will be improved (might 
30 % efficiency). Then, among 10 s BB events, ~ 10 5 events without any constraint 
on nix may be reconstructed. For more details on the problems of experimental 
reconstruction and continuum background, see Ref. ||. 

Even without a full event reconstruction, one may have good measurements of 
missing energy and momentum of missing neutrino, E v and p v satisfying the mass- 
shell condition, 

E 2 u -f v = ml = Q. 

In this case, the hadronic mass nix is not fully constructed, but it is bounded by 

m x < m x,max =m 2 B + mf u - 2^m B (E lu - (3p lv ). (20) 

Here, m B is the B— meson mass, and 7 = (1 — /5 2 ) _1//2 = my/2m B for the symmetric 
B factory with e+e" -> T(4£) -> BB. Since (3 is very small, m x max is close to m 2 x , 
and we lose very little efficiency It turns out that ~ 80% events for b —>■ u transitions 
survive the cut on the m 2 x max : 

m x .mux 

< m D . 

In any case, studying the hadronic mass distributions in inclusive semileptonic 
b — » u transition is experimentally viable. It's also theoretically better described, 
so theoretical uncertainties in determining \V u b\ would be less compared to the \V u b\ 
determined from studies of exclusive decay modes. In summary, we would have better 
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statistics in extracting \V u b\ by measuring the mx distributions in inclusive b — > u 
semileptonic decays, and have better theoretical handles over the inclusive decays 
rather than exclusive decays. 
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Figure Captions 

Fig. 1 The m x distributions in B — > X CjU lu with \V ub /V cb \ = 1. The b — > c 
transition is dominated by the X c = D, D*, D**, and can be reliably calculated in the 
HQET. The regions between the arrows are predicted rate in the unit of 10 13 sec -1 
when the dmx integration over the delta function is performed. On the other hand, 
the b — > u transition is largely nonresonant. The cases with X u = ir,p are shown 
explicitly for two different models. For X u = 7r, the region between the open triangles 
are predictions by Hagiwara et ai, and the region between the closed triangles are 
predictions by Casalbuoni et ai. For X u = p, the regions between the lower two 
curves are predictions by Hagiwara et ai, whereas the region between the upper two 
curves are predictions by Casalbuoni et ai. The inclusive m x distribution for b — > u 
was obtained from the ACCMM model with hadronic mass constraint of mx ~ Ivn^. 
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